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I. INTRODUCTION 



Molecular vibrational and rotational spectroscopy is one of the important 
parts of molecular physics and one of the main tools for other scientific areas 
such as biology [1] and environmental sciences [2]. The harmonic oscillator 
could be useful ground to explain the molecular vibrations but this model 
is restricted for only Lowest states [3]. To improve the theory of molecular 
vibrations the anharmonic oscillators such as Morse and Mie-type potentials 
(Kratzer potential and its generalization) can be used to solve exactly the 
Schrodinger equation (SE) and provide more reliable model for diatomic 
molecules [4]. Mie-type potentials being considered as an example in the 
present work have important advantages such as having eigenfunctions behav- 
ing correctly at r — > and r — > oo and providing exact solutions to the SE 
[5]. So, these potentials have been used to determine molecular structures 
and received much attention in Literature [6]. In the present work, we deal 
with also another diatomic potential called pseudoharmonic potential pro- 
posed by Davidson [7]. This potential is used to describe the roto-vibrational 
states of diatomic molecules and nuclear rotations and vibrations [8]. 

In the light of the above considerations, it could be interesting to solve 
exactly the SE for the pseudoharmonic potential and Mie-type potentials 
and find any Estate solutions in the view of molecular physics phenomenon. 
Moreover, obtaining the exact solutions of the SE for the molecular poten- 
tials is one of the main problems in quantum physics [6]. One of the methods 
giving exact solutions of the SE is used in Ref. [9] where the wave equation 
is solved for the non-central potential within the framework of the supersym- 
metric quantum mechanics. In Ref. [10], the energy eigenvalues of the radial 
SE are obtained for the Coulomb potential by using path integral formal- 
ism and the author also stated how can be obtained the wave functions. In 
Ref. [11], energy spectrum of the Coulomb, Morse and harmonic oscillator 
potentials have been studied by using point canonical transformation where 
the formalism has been extended to the case of position-dependent mass. 

We List some methods used in Literature to solve the wave equations for the 
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pseudoharmonic potential and Mie-type potentials: Nikiforov-Uvarov method 
[12, 13], algebraic approach [14], polynomial solution [15], exact quanti- 
zation rule [16], hypervirial theorem with perturbation theory [17], shape- 
invariance procedure [18], solutions in terms of hypergeometric functions 
[19], etc.. In this work, we find exact bound state solutions of pseudohar- 
monic potential and Mie-type potentials by reducing the SE to a first-order 
differential equation via Laplace transform approach (LTA) and therefore we 
make use of integral to the obtain energy eigenvalues and the corresponding 
eigenfunctions. Actually, the LTA is an integral transform which has been 
used by many authors to solve the SE for different potentials [20, 21, 22, 
23]. The LTA could be a nearly new formalism in the Literature and serve 
as a powerful algebraic treatment for solving the second-order differential 
equations. As a result, the LT methods describe a simple way for solving 
of radial and one-dimensional differential equations. The other advantage 
of this approach is that a second-order equation can be converted into more 
simpler form whose solutions may be obtained easily [22]. 



II. ENERGY EIGENVALUE SOLUTIONS 

Time-independent Schrodinger equation is written as 

j- J^-V 2 + V(r) j vi/(r, 0, cp) = E ni W(r, 6, ip) , (1) 

and defining the wave function M/(r, 6, ip) — ^R(r)Y(0, cp), we obtain the 
radial SE as [6] 

{£-^ + £[*.-™]}KCr)-0. 

where I is the angular momentum quantum number, m is the particle mass 
moving in the potential field V(r) and E nl is the nonrelativistic energy of 
particle. 
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A. Pseudoharmonic Potential 



The pseudoharmonic potential is given [8] 

V(r) = cur 2 + + a 3 , (3) 
r 2 

where a*(i = 1,2,3) are real parameters. Inserting Eq. (3) into Eq. (2), we 
obtain 



where 

2 2raai . . 2ma 2 , „,„ , , ^ 2 2m . , , 

M 2 ; u(u + 1) = + ^(^ + 1) ; £ 2 = — (E n , - a 3 ) . (5) 



Defining the new variable y — r 2 and rewriting the radial wave function as 
R(y) — y- l//2 0(y), Eq. (4) turns into 

By using the Laplace transform defined as [24] 

£ {0(V)} = fit) = / dye- ty 0(y) , (7) 

Jo 

Eq. (6) reads 

which is a first-order ordinary differential equation and its solution is simply 
given 

/ »\ -£-*("+§)/ »\ £-*("+§) 
/(*) = A/(* + f) ■ O) 

where A/ is a integral constant. In order to obtain finite wave functions, it 
should be 



— " 2 + 2 ' = U ' ( n = ' 1 ' 2 ' 3 ----) ( 1C} ) 



4 



which gives single-valued wave functions. By using this requirement and 
expanding Eq. (9) into series, we get 

" (-l)*n! (t + H) 

f(t) = N'\ - — —— - — y- , en) 

fc=0 

where N' is a constant. By using the inverse Laplace transformation [24] we 
immediately obtain the solution of Eq. (6) 

0(y) = n" f ( ~ 1)fcn! r( ^| f) y(^) e -^ , (12) 



fc=0 

/// 



where N" is a constant. On the other hand, the confluent hypergeometric 
functions is defined as a series expansion [25] 



(n — m)!m! T(cr + m) 

m=0 

So, comparing Eq. (12) with Eq. (13) we deduce that 



0(y) = A/"'e-'«"V + * iFi(-n, i/ + |, i/) , (14) 
We obtain finally the radial wave functions 

Rty) = j^ e -^/2y(u+i)/2 lFl( _ n> t/ + |, t/) . (15) 

where A/" is normalization constant. Using the normalization condition given 
as / °° [R(r)] 2 dr — 1 and the relation between the Laguerre polynomials and 
confluent hypergeometric functions as L^{x) — ^irtp+i? lFl (~ n < V + 
[25] , the normalization constant in Eq. (15) is written 



M = M ("+3/2)/2. 

where we have used [25] 



f 2r(n + ^ + f) 
n! 



(16) 



/" x«e- x L q n (x)L q n ,(x)dx = F(Q + ™ + 1} fi nn , . (17) 

In Fig. (1) we show the variation of the normalized wave functions of 
the pseudoharmonic potential on the coordinate r. We give first six wave 



functions according to the quantum number pairs (n, I). Inserting the param- 
eters in Eq. (5) into Eq. (10), the energy spectrum of the pseudoharmonic 
potential is obtained 

£ "' = ^ + \/^(" + l + W 1+4£(£ + 1) + ^) ■ (18) 

We give our numerical energy eigenvalues for two different diatomic poten- 
tials in Table I. We compare our results with the ones given in Ref. [12] by 
setting the potential parameters as ai = a 2 = D r 2 and a 3 = — 2D . 
Let us study the results of some special cases. Firstly, if we put a 3 = we 
obtain from Eq. (18) 

,_ /8ft 2 ai / ,1,1 [~ , nfn , ~7 , 8ma 2 \ 

E "' = V— ^ n+ 2 + 3V 1+4£(£+1) + ^J ' (19) 

which is the same result given in Ref. [26] for the potential of the form 
£ + Br 2 . Secondly, if we choose the parameters as a 2 = a 3 = and 
ax — \mijj 2 in Eq. (18) we get 

Ent = TiuJ ^2n + I + |^ , (20) 
where if we define n' — 2n + £ as 'principal quantum number' we obtain 

E n , t = Tiu (n' + |^ . (21) 
which is exactly the spectrum of three-dimensional harmonic oscillator [26]. 

B. Mie-type Potentials 

The Morse potential is an example of this type of potentials or Kratzer 
potential and its generalization having the forms, respectively, 

V(r) = -0 - if) , (22) 

and 

V(r) = D(^—^y . (23) 
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where D is the dissociation energy and r is the equilibrium distance [15]. 
So, the Mie-type potentials can be simply given as 

W> = £ + £ + C, (24) 

where a,b,c are real potential parameters. Inserting Eq. (24) into Eq. (2), 
redefining the wave function as R(r) — Vrip{r) and using the following 
abbreviations 

7 2 — — + ^ + 1) + - ; <5 2 = — — ; e 2 — (c - E ni ) , (25) 
Ti 4 71 Ti 

gives 

( r2 ^r + r ^ - [4 + - + £2 1 r2 l ^ = ■ < 26 > 

t dr 2 dr ir 2 r J J 
Setting </?(r) = r a 0(r) with a is a constant and then inserting into Eq. (26) 
Leads 

{ r 2 ^- + (2a + l)r— - £ 2 r 2 - <5 2 r + a 2 - -y 2 l 0(7") = , (27) 
t dr 2 dr J 

In order to obtain a finite wave function when r — > oo, we must take a = —7 
in Eq. (27) and then we get 

- (27 - 1)^ - <5 2 - £ 2 r| 0(r) = . (28) 

Applying the Laplace transform to Eq. (28) we obtain a first-order differen- 
tial equation 

(t 2 - E 2 ) ^ + [(27 + 1) t + <5 2 ] /(t) = , (29) 
whose solution is 

_ 2-1+1 

f(t) = N (t + £) - (27+1) (1^1) * 2 ■ (30) 

The wave functions must be single-valued which requiring that 

<5 2 27 4- 1 

- 4?— = n, (n = 0., 1,2,3, .. .) (31) 

2e 2 
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Taking into account this requirement and applying a simple series expansion 
to Eq. (30) gives 



n 



fc=0 



(32) 



where /V 7 is a constant. Using the inverse Laplace transformation [24] in Eq. 
(32) we deduce that 



0(r) = N"r^e- er Y -± — 

<• (n — 

Finally, we obtain 



(-l) fc n! r(27+l) 



fc=0 



(n - fc)!/c! T(27 + 1 + /c) 



(2er) fe , 



(33) 



<p(r) = T (-D fc n! r(27 + D (2£r)fc 

^ v y (n - /c)!/c! f(27 + 1 + /c) 



(34) 



fc=0 



Comparing last equation with Eq. (13) we write the radial wave function as 

R(r) = Mr^e- £r 1 F 1 (-n, 2-y + 1, 2er) . (35) 

where the normalization constant is given by following the same procedure 
in previous section as 



M= r(27+ !)• 



'n!(2n + 27 + I) 



(36) 



r(n + 27 + 1) 

We give the dependence of the wave functions of the Kratzer potential on 
spatially coordinate r in Fig. (2) where the wave functions are plotted for the 
same quantum number values as in pseudoharmonic potential. Using Eqs. 
(25) and (31) we obtain the energy eigenvalues of the Mie-type potentials 



Ene — c — 



8m 



-i 2 



2mb/h 7 



n + \ (l + 2> /2=£ + *(*+!) + J ) 



(37) 



which is the same result with the ones obtained in Ref. [27]. We summarize 
our numerical results for different quantum number pairs (n,£) in Table II. 
To compare our results we chose the potential parameters as a — D e r\, b — 
—2D e r e , c—D e used in Ref. [13]. 
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III. RESULTS 



Our numerical energy eigenvalues of two diatomic molecules interacting 
in short-range [28] given in Tables I and II have a good accuracy with the 
ones obtained in Refs. [15] and [13]. Figs. 1 and 2 show the variation of 
the wave functions versus r for the pseudoharmonic and Kratzer potentials. 
The wave functions of the pseudoharmonic (Kratzer potential) go to zero as 
r — y and as r ~ 8 (r ~ 35). The eigenfunctions corresponding to n — 1 
for each potential go to zero from lower part of the vertical axes while the 
remaining functions from upper part of the zero axes. 

IV. CONCLUSIONS 

We have exactly solved the radial Schrodinger equation for the pseudo- 
harmonic and Mie-type potentials by using Laplace transform approach. We 
have found the energy eigenvalues and the corresponding normalized eigen- 
functions of the diatomic potentials. We discussed briefly some special cases 
of the potentials. We observed that our analytical results and also the re- 
sults for the special cases are the same with the ones obtained in literature. 
We also summarized our numerical energy eigenvalues for two different di- 
atomic molecules. It seems that the Laplace transform approach is very 
economical method about solving the wave equations for some potentials. 
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TABLE I: Energy eigenvalues of the pseudoharmonic potential for different 
values of n and I in eV (the parameter values are used in Ref. [12]: 

D =96288.03528 cm -1 , r Q =1.0940 A, m =7.00335 amu for A/ 2 ; 
D =87471.42567 cm" 1 , r =1.1282 A, m =6.860586 amu for CO). 



A/ 2 CO 



n 


I 


Our ResuLts 


Ref. [15] 


Our Results 


Ref. [15] 








0. 


,109180 


0.1091559 





,101953 


0.1019306 


1 





0. 


,327414 


0.3273430 





,305738 


0.3056722 




1 


0. 


.327913 


0.3278417 





.306217 


0.3061508 


2 





0. 


,545648 


0.5455302 





.509524 


0.5094137 




1 


0. 


,546147 


0.5460288 





.510003 


0.5098923 




2 


0. 


,547145 


0.5470260 





.510961 


0.5108495 


3 





0. 


,763883 







,713310 






1 


0. 


,764382 







,713789 






2 


0. 


,765380 







,714747 






3 


0. 


,766877 







.716183 




4 





0. 


.982117 


0.9819045 





.917095 


0.9168969 




1 


0. 


.982616 


0.9824031 





.917574 


0.9173755 




2 


0. 


,983614 


0.9834003 





.918532 


0.9183327 




3 





.985111 


0.9848961 





,919969 


0.9197684 




4 


0. 


,987107 


0.9868903 





,921885 


0.9216825 
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TABLE II: Energy eigenvalues of the Kratzer potential for different values 

of n and I in eV . 



N 2 CO 



n 


I 


Our Results 


Ref. [13] 


Our Results 


Ref. [13] 








0.054434 


0.054430 


0.050827 


0.050823 


1 





0.162068 


0.162057 


0.151296 


0.151287 




1 


0.162557 


0.162546 


0.151765 


0.151755 


2 





0.268245 


0.268229 


0.250369 


0.250354 




1 


0.268728 


0.268711 


0.250831 


0.250816 




2 


0.269692 


0.269675 


0.251756 


0.251744 


3 





0.372992 


0.372972 


0.348070 


0.348051 




1 


r\ 070 /i /c o 
U.o f o4oo 


(J.o f o44f 


r\ o a o cz 
U.o4o52o 


U.o4ob(J ( 




2 


0.374419 


0.374398 


0.349438 


0.349418 




3 


0.375846 


0.375823 


0.350806 


0.350785 


4 





0.476334 


0.476313 


0.444425 


0.444403 




1 


0.476803 


0.476779 


0.444871 


0.444852 




2 


0.477742 


0.477717 


0.445774 


0.445751 




3 


0.479150 


0.479124 


0.447123 


0.447099 




4 


0.481026 


0.480999 


0.448921 


0.448895 
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IG. 1: Variation of the first six normalized eigenfunctions of the 
pseudoharmonic potential (in m = h — 1 unit). 
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FIG. 2: Dependence of the first six normalized eigenfunctions of the 
Kratzer potential on spatially coordinate r (in m = ft = 1 unit). 
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